First principles electronic structure of spinel LiCr204: A possible half-metal ? 



O 

o 



in 



C/2 



O 

(N 
> 

O 
cn 
o 



I 

O 

o 



Markus LaueiQ 
Fachrichtung Theoretische Physik, Universitat des Saarlandes, 
Postfach 15 11 50, D-66O4I Saarbriicken, Germany 

Roser ValentQ 
Institut fiir Theoretische Physik, Universitat Frankfurt, 
Robert-Mayer-Strasse 8, D-60054 Frankfurt, Germany 

H. C. Kandpai 

Institut fiir Anorganische Ghemie, Universitat Mainz, Duesberg Weg IO-I4, D-55099 Mainz, Germany 

Ram SeshadrJl 

Materials Department, University of Galifornia, Santa Barbara, CA 93106 USA 

We have employed first-principles electronic structure calculations to examine the hypothetical 
(but plausible) oxide spinel, LiCr2 04 with the d^'^ electronic configuration. The cell (cubic) and 
internal (oxygen position) structural parameters have been obtained for this compound through 
structural relaxation in the first-principles framework. Within the one-electron band picture, we 
find that LiCr2 04 is magnetic, and a candidate half- metal. The electronic structure is substantially 
different from the closely related and well known rutile half- metal Cr02. In particular, we find a 
smaller conduction band width in the spinel compound, perhaps as a result of the distinct topology 
of the spinel crystal structure, and the reduced oxidation state. The magnetism and half-metallicity 
of LiCr204 has been mapped in the parameter space of its cubic crystal structure. Comparisons 
with superconducting LiTi204 (d'^"^), heavy-fermion LiV204 (d^'^) and charge-ordering LiMn204 
{d^'^) suggest the effectiveness of a nearly-rigid band picture involving simple shifts of the position 
of Ey in these very different materials. Comparisons are also made with the electronic structure 
of ZnV204 (d^), a correlated insulator that undergoes a structural and antiferromagnetic phase 
transition. 



PACS numbers: 71.15.Nc, 71.20.-b, 75.50.- 



I. INTRODUCTION 

The stoichiometric, self-doped oxide spinels LiM2 04 
(with M = Ti, V and Mn) can be said to constitute 
a periodic "hall of fame". In all of them, Li is mono- 
valent and the oxygen ion, negative and divalent; this 
means that the transition metal oxidation state is for- 
mally 3.5+. The number of d electrons per transition 
metal atom would therefore be 0.5, 1.5, 2.5 or 3.5 when 
M = Ti, V, Cr or Mn respectively. LiTi204 is a super- 
conductor with a Tc of 12 K. When first prepared in 1973 
it was described as a high temperature superconduc- 
tor. Along with BaBii_a;Pba;03 2], it was amongst the 
only oxide superconductors with Tq > 10 K until the ad- 
vent of the layered cuprates. It is also the only metallic 
oxide spinel known. Its neighbor, LiV204 is already on 
the verge of localization, and as a result, displays heavy- 
fermion behavior j^] — one of the few systems without / 
electrons that is known to do so — LiMn2 04, one more 
member from this family, is not only a very important 
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cathode material in secondary Li-ion batteries ^ , but 
also displays charge-ordering [3, , associated with half- 
integral charge that can order over two distinct sites. In- 
deed, the well-known Verweij transition in ferrite associ- 
ated with the ordering of charge on the spinel lattice 
can perhaps be best understood by examining LiMn204. 
Near the charge-ordering temperature of this compound, 
a magnetic field can influence electrical transport and 
this translates to negative magnetoresistence 0- 

Conspicuous in this list by its absence is the spinel 
LiCr204. Attempts to prepare this compound are 
fraught with the difficulty of stabilizing Cr in a formal 
oxidation state that is somewhere between 3-f and 4-1- . 
Indeed ferromagnetic Cr02 with the rutile structure is 
well known to require high oxygen partial pressures for 
its preparation. Cr02 decomposes at ambient oxygen 
partial pressures when the temperature exceeds 473 K 
Given that (i) the isostructural compounds neigh- 
boring LiCr204 are so interesting, and (ii) the structure 
and oxidation state of LiCr204 are reminiscent of Cr02, 
the prototypic oxide half-metal 0, 0| , we thought the 
compound LiCr204 worthy of study. 

In this work, we use spin-polarized first principles den- 
sity functional calculations within two different schemes: 
Full potential linearized augmented plane-wave calcula- 
tions [T^ [1^ that permit an optimization of the geome- 



try of the cubic spinel compound (lattice parameter and 
internal structural parameter) and Linear Muffin Tin Or- 
bital (LMTO) calculations which are geared to the 
visualization of spin-resolved bonding. The electronic 
structure of LiCr204 is compared with that of rutile 
Cr02, LiM204 (M = Ti, V and Mn), and ZnV204. The 
parameter space of half-metallicity in LiCr2 04 is mapped 
out in terms of the two structural parameters. 



II. CRYSTAL STRUCTURES AND 
COMPUTATIONAL METHODS 

The cubic spinel crystal structure AB2O4 is completely 
described by two parameters, the cubic cell parameter 
a, and the internal oxygen positional parameter 
where x ~ 0.25 The A atoms (Li) occupy tetrahe- 
dral sites created by a network of BOe octahedra where 
B is here a transition metal. The octahedra of oxygen 
around B are perfectly regular when x = 0.25. Devia- 
tions 5, X = 0.25 -|- 5 indicate either a trigonal expansion 
{6 > 0) or compression {S < 0), achieved by pushing 
two opposite triangular faces of the octahedron apart or 
together. The B atom sublattice (Cr in LiCr204) com- 
prises a network of B4 tetrahedra. Because the B atoms 
form tetrahedra, antiferromagnetic ordering is frustrated 
on the B site of the spinel 16] in a manner reminiscent of 
the residual entropy problem in the crystal structure of 
ice Ih |r^. Various views of the spinel crystal structure, 
and the structure of rutile Cr02 are shown in Fig. ^ 

Scalar-relativistic Kohn-Sham equations were solved 
taking all relativistic effects into account except for the 
spin-orbit coupling. We performed both full potential 
linearized augmented plane-wave (LAPW) calculations 
based on the WIEN97 and WIEN2K codes [H El 
as well as Linear muffin-tin orbital (LMTO) calcula- 
tions jl4j within the atomic sphere approximation (ASA) 
where we used the Stuttgart tb-lmto-asa program 
[rsf. The generalized gradient approximation (GGA) for 
the calculation of exchange correlation was considered 
following the Perdew-Burke-Ernzerhof |l9j | parameteri- 
zation for the LAPW calculations and the Perdew-Wang 
prescription 20] for the LMTO calculation[2H. 

In our LAPW calculations we considered twenty in- 
equivalent sampling k points and the modified tetrahe- 
dron method 22] for the Brillouin zone integration during 
the self-consistent iterations. In order to test the accu- 
racy of the sampling, we also performed calculations with 
up to 120 and 195 irreducible fc-points without observing 
significant qualitative changes in the results apart that 
some spikes in the density of states plots get rounded. 
We set the energy threshold between core and valence 
states at -6 Ryd 23J. For the number of plane waves, 
the criterion used was RMT (Muffin-Tin Radius) x fcmax 
(for the plane waves) = 8. We considered various sets of 
muffin-tin radii to ensure well converged spin-polarized 
calculations. 

In our LMTO calculations 256 irreducible k points 
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FIG. 1; (a) Spinel AB2O4 structure showing octahedral B 
atoms (grey spheres) at the centers of Oe octahedra (depicted 
using sticks). The black spheres are the tetrahedrally coordi- 
nated A atoms, (b) Four BOe octahedra in the spinel struc- 
ture showing the formation of B4 tetrahedra as well as the 
nature of the edge-sharing of oxygen, (c) Network of corner- 
sharing B4 tetrahedra in the spinel structure, (d) Portion of 
the crystal structure of rutile Cr02 showing Cr atoms (grey 
spheres) in the middle of Oe octahedra. 



were used in the primitive wedge of the BZ. Space filling 
in the unit cell of LiCr204 was ensured through the use 
of two empty spheres with bases of s and p orbitals at (0, 
0, 0) and {y, 1/8, 1/8) where the value of y is set auto- 
matically by the code according to the oxygen position. 
LMTO electronic structures were analyzed by calcula ting 
crystal orbital Hamiltonian populations (COHPs) |2J] 
which are densities of states weighted by appropriate 
Hamiltonian matrix elements. COHPs are indicative of 
the strength and nature of a bonding (positive COHP) 
or antibonding (negative COHP) interaction. The signs 
we use here are the opposite of what is used in the orig- 
inal definition of Dronskowski and Blochl 0]. LMTO- 
ASA calculations do not allow for structural relaxation 
while structural relaxation process is implemented in the 
LAPW package. We will present in what follows a com- 
parison of the results obtained by both methods. 
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III. CRYSTAL AND ELECTRONIC 
STRUCTURE OF LiCr204 

A. Structure optimization 

Within the LAPW scheme we performed a structure 
optimization for LiCr204 where the cubic cell parameter 
a and the oxygen position x were varied till the optimal 
structure was obtained. In order to perform the opti- 
mization we considered two steps: (i) We first took sev- 
eral crystal structures with different lattice constants a 
and the same x. For each of them we performed an iter- 
ation procedure up to self-consistency and we compared 
the total energies. The energy minimum then defines the 
optimal a. For LiCr204 the energy minimum was reached 
at a = 8.11(4) A (see Fig.|2). 
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FIG. 2: LAPW energy versus primitive unit cell volume (in 
Bohr"^) for the LiCr2 04 system. The squares are calculated 
LAPW data and the dashed curve is a fit to the Murnaghan 
equation of state 

(ii) In a second step we proceeded to optimize the 
oxygen positions by evaluating the atomic forces [26j |. 
The equilibrium oxygen position was obtained using a 
damped Newton dynamics method according to the ex- 
pression 
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where i?^ and are the coordinate and force for the 
atom m at time step t. 5m determines the speed of motion 
and r^TO changes from rjm to 1 — rjm if the force changes its 
direction from one step to the next. The optimal value we 
obtained for the oxygen position was x= 0.253(5) which 
implies a small distortion from the perfect spinel struc- 
ture. A further unit-cell optimization did not change the 
result significantly. 

In the previous calculation scheme the unit-cell is first 
optimized and then the structure relaxed. Since the two 
operations do not necessarily commute, we considered 
also first a structural relaxation in a bigger unit cell fol- 
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FIG. 3: (a) LAPW and (b) LMTO projected densities of 
states for LiCr2 04. Black traces indicate Cr d states and 
grey traces indicate O p states. Spin up states are shown in 
the upper portions and spin down states in the lower portions 
in both panels. The origin on the energy axis in this and in 
the following plots is the E-p- 



lowed by a lattice constant optimization. The two proce- 
dures do not yield exactly the same optimal value for a 
and X, but deviations from a = 8.11(4) and x — 0.253(5) 
are small and can be considered within the error bars. 

We checked both in LAPW and LMTO that the ferro- 
magnetic spin-polarized configuration is lower in energy 
than the non-spin-polarizcd case. We also considered 
within LAPW an antifcrromagnetic spin arrangement of 
the Cr atoms in the spinel lattice which showed to be 
energetically unfavourable with respect to the ferromag- 
netic arrangement. 



B. Density of states (DOS) 

Figure |31 compares the partial density of chromium d 
states and oxygen p states from spin polarized (a) LAPW 
calculations and (b) LMTO calculations. In both plots, 
the Fermi energy is taken as the origin on the energy 
axis. The upper panel of each plot displays majority 
(up) spin states and the lower panel minority (down) 
spin states. Cr-d states are shown using black traces and 
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0-p state using grey traces. The optimized crystal struc- 
ture obtained from the LAPW calculations was used as 
the input crystal structure for the (unrelaxed) LMTO 
calculations. We note immediately that in terms of gross 
features, the two computational methods yield quite sim- 
ilar results, in terms both of band positions and widths, 
as well as details in the structures of the projected den- 
sities of state. Differences in the two plots arise only as 
a result of slightly larger gaps between occupied and un- 
occupied states between the two computational schemes; 
the LAPW calculations indicate larger (by 0.5 eV) crys- 
tal field splitting as well as larger exchange splitting (by 
the same magnitude). This different energy shift can be 
understood in terms of the different nature of the two 
schemes and that for open structures -as the spinels are- 
the LMTO approach is not as accurate as the LAPW 
approach. We use here therefore the LAPW results as a 
reference for accuracy. 

O p states are found to be largely centered around -4.5 
eV and possess a width of about 5 eV. There are also 
Cr d states in this energy region, indicative of covalency. 
Near Ep we find majority Cr t2g (t) and no minority 
spins in the LAPW DOS, i.e. the polarization is 100% 
what characterizes the system as half-metallic. In the 
LMTO DOS there is a residual spin density of minor- 
ity Cr t2g (i) spins at Ep. While strictly speaking the 
LMTO results predict an approximate half-metal behav- 
ior, the small overall differences between both LMTO 
and LAPW approaches justifies the use of the LMTO- 
ASA calculations in some of the subsequent discussion. 

The t2g{1) states have a width of about 2 eV and are 
bimodal, with the two modes separated by a pseudo-gap. 
The Ep lies in the center of a peak in the DOS. The crys- 
tal field splitting between t2g(T) and eg(t) is of the order 
of 3 eV (taking the midpoints of the different manifolds) 
and the exchange splitting between i2g(T) and i2g(J.) is 
of the order of 2 eV. These values are in the LMTO DOS 
-as discussed above- about 0.5 eV larger. 

The two-peaked t2g(T) DOS is puzzling. One suspect is 
the trigonal distortion of the oxygen octahedron around 
the Cr atoms. In LiCr204, the calculated oxygen posi- 
tion 0.253(5) suggests a very small distortion from perfect 
octahedral coordination. We have performed LMTO cal- 
culations using the ideal oxygen position x — 0.250 cor- 
responding to perfect CrOe octahedra. The differences 
in the DOS (between the structures where x = 0.25 and 
X = 0.253) are negligible and the two-peaked structure 
of the i2g(T) DOS is retained, and Ep continues to fall 
in the second DOS peak even when O atoms are in the 
ideal position and the CrOg octahedra are undistorted. 

Bimodal t2g(T) states, also observed in the electronic 
structure of spinel LiTi2 04 27], could arise from the re- 
stricted B-O-B dispersion in the spinel structure; this 
network can be built up from B4O4 cubes that are in- 
terconnected through B atoms. The B-O-B bond angles 
are therefore 90°. This limits dispersion significantly. In 
fact, through the analysis of the COHP in LiCr204 in the 
next section we can identify that the upper lobe of the 



DOS corresponds mainly to Cr-0 interactions and less to 
Cr-Cr interactions. 

The most important conclusion that we draw from the 
DOS of ferromagnetic LiCr204 is that states at and near 
E-p are strongly spin-polarized and the system is indeed a 
magnetic half-metal. However, the fact that Ep falls on 
a sharply peaked region of the DOS suggests (i) that the 
compound might be unstable/difficult to prepare and (ii) 
if it is prepared, might be subject to electronic instabili- 
ties associated with electron correlation, i.e. the opening 
of a Mott-Hubbard gap, or charge-ordering. We suggest 
that various different scenarios could stabilize the struc- 
ture of LiCr204, for instance the presence of a slightly 
distorted antiferromagnetic structure as it is the case in 
ZnV204. Also if the system gets a little oxidized, by 
removing for instance some Li to form Li5Cr2 04, where 
(5 < 1, it may be then possible to have Ep sitting in 
the pseudogap of the DOS, i.e. at the center of the bi- 
modal DOS and therefore to have a stable structure. A 
third possibility, as mentioned above, is the explicit con- 
sideration of the electron correlation which may open a 
Mott-Hubbard gap at Ep and therefore settle stability. 



C. Comparisons with Cr02 

It is important to compare the electronic structure 
of LiCr204 with the prototypic oxide half metal Cr02. 
Present understanding of this material has been reviewed 
recently This rutile compound Fig. dd)) has (P 

Cr^+ surrounded octahedrally by O^^. The octahedra 
form ribbons that share edges, and each Cr can therefore 
be said to bond with two others through the shared edge 
with a Cr-Cr distance of 2.92 A A number of au- 

thors 30] have described the electronic structure of this 
material, including Korotin et al. jsil ] who performed 
LDA + U calculations and argue thereof that the one- 
electron description of a ferromagnetic band metal does 
not suffice and that U — 3.0 eV is required to explain the 
essential features of this compound. Mazin et al. [s^ ar- 
gued that when the local interaction is small compared 
to the bandwidth (W) U/W« 0.5 as in Cr02, it is dis- 
putable which method spin-polarized LDA or LDA+U 
provides a better description of the system. 

We have performed spin-polarized LDA calculations in 
the LMTO scheme on Cr02 for two reasons. The first is 
that it provides us a consistent comparison with LiCr204. 
The second is that spin-polarized bonding through the 
use of COHPs has not been presented before for this im- 
portant compound. Figure ^compares the Cr d and O p 
DOS of (a) LiCr2 04 and (b) Cr02. O p states of CrOa 
are found to be quite similar to what is seen for LiCr204, 
spreading from around -7 eV to -2 eV. Once again, there 
is a strong admixture of Cr d states in this energy region. 
The region near the Fermi energy comprised largely of 
Cr d states. The crystal field splitting between i2g(T) 
and eg{]) is also similar in both compounds, as might 
be expected considering the similarity in charge and co- 
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FIG. 4: LMTO projected densities of states for (a) LiCr204 
and (b) Cr02. Black traces indicate Cr d states and grey 
traces indicate O p states. 



ordination. What is significantly different is the larger 
t2g{1) bandwidth in Cr02, as a result of the different Cr- 
O topology in the rutile structure, that allows for eas- 
ier metal-oxygen-metal hopping. Unlike what we see in 
LiCr204, the in Cr02 does not lay on a peak but 
instead, hes on the edge of a pseudo-gap in the i2g(T) 
states. This gap arises due to a two-in, fqur-out Jahn- 
Teller distortion of the CrOg octahedra [23 . 

Crystal orbital Hamiltonian populations (COHPs) per- 
mit a better understanding of the precise nature of the 
different states that would determine the properties of 
LiCr204. In Fig. we plot (a) Cr-Cr COHPs and (a') 
Cr-0 COHPs as a function of energy for LiCr2 04, per 
primitive unit cell. This means Cr interacts with six O 
(at a distance of 2.00 A) and to six other Cr (at a dis- 
tance of 2.86 A). The COHPs are spin resolved into inter- 
actions between orbitals with majority spin, and orbitals 
with minority spin. In the absence of spin-orbit coupling, 
majority states cannot bond with minority states. In the 
region of the DOS which was predominantly 0-p (-7 eV to 
-2 eV) we find evidence for quite strong Cr-0 covalency. 
At the E-p, we find a combination of majority antibonding 
Cr-Cr interactions and majority antibonding Cr-0 inter- 
actions, being the Cr-0 interaction predominant. There 
are no minority interactions near the E-p suggesting that 



the bonding is highly spin-polarized. This is what makes 
this compound a putative half-metal. 

Similar COHPs for Cr02 [Fig. E^b) and (b')] are de- 
picted per unit cell. In Cr02, each Cr has 6 oxygen 
neighbors (4 at 1.90 A and 2 at 1.91 A), but only two 
Cr neighbors (at 2.92 A). The shorter Cr-0 distance in 
Cr02, compared with LiCr204 reflects the slightly higher 
oxidation state of Cr. The Cr-0 COHP suggests strong 
Cr-0 bonding in both spin directions in the region where 
O p states are found. Cr-Cr bonding is insignificant in 
this structure, as a result of the larger Cr-Cr distance 
than what is found in the spinel, as well as fewer neigh- 
bors. This means that Cr t2g{'\) states at the Ep are 
largely non-bonding. As in LiCr204, the Cr-0 antibond- 
ing states at Ep are completely spin-polarized and cor- 
respond to interactions in the majority spin channel. 



D. Comparisons with other LiM204, M = Ti, V 
and Mn, and with ZnV204 

It is of interest to compare the DOS of the compounds 
LiM204, M = Ti, V, Cr and Mn. LiTi204 is a super- 
conductor. LiV2 04 is a metal on the verge of a metal- 
insulator transition, with correspondingly, an unusually 
large effective carrier mass LiMn204 is a correlated 
antiferromagnetic insulator at low temperatures and dis- 
plays a charge-ordering transition near the room tem- 
perature 0, la 111 • A number of recent papers describe 
the electronic structure of LiV204 l33l|. For example, 
Singh et al. ,34] have performed LAPW calculations on 
LiV204. They find some evidence for a separation of the 
t2g manifold into aig and states as a result of the trigo- 
nal distortion of the VOe octahedra, which is large in this 
compound. This results in the presence of both flat and 
disperse bands near the Ep^ the flat bands presumably 
giving rise to the effective mass enhancement. Mishra 
and Ceder "sS] have reported density functional calcula- 
tions on LiMn204, but with an emphasis on structural 
stability, rather than on magnetism. 

Figure compares the LMTO DOS of LiM204, M 
= Ti, V, Cr, Mn. For the V, Cr and Mn compounds, 
we found that the ferromagnetic spin-polarized configu- 
ration was energetically more favourable that the non- 
spin-polarized configuration. For the Ti compound, the 
ferromagnetic solution lies in LAPW and LMTO higher 
in energy than the non-spin-polarized solution. In fact in 
LMTO we started for this system with a ferromagnetic 
spin-polarized configuration and in the self-consistency 
cycles it gets increasingly non-spin-polarized. The width 
of the O p states in all three compounds is similar. The 
separation between O p and metal d states is largest in 
the V compound and smallest in the Mn compound. This 
is in keeping with the expectation that for a given oxida- 
tion state and coordination, moving to the right amongst 
transition metals results in a stabilization of metal d 
states. Such stabilization is in fact the basis of the evolu- 
tion of the band gap in correlated transition metal com- 
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FIG. 5: Crystal orbital Hamiltonian populations (COHPs) in LiCr2 04 and Cr02 per primitive unit cell, (a) Cr-Cr interactions 
in LiCr204, (a') Cr-0 interactions in LiCr204,(b) Cr-Cr interactions in Cr02, (b') Cr-0 interactions in Cr02. In each panel, 
the black trace corresponds to interactions in the spin up channel and the grey trace to interactions in the spin down channel. 



pounds; from the Hubbard U type (arising from d — d 
Coulomb correlation) for the early transition metals to 
charge-transfer A (arising from ligand to metal charge 
transfer) in the later transition metals [s^ls^l . The width 
and gross features of metal t2g states in these compounds 
display many similarities, and in fact, a simple rigid band 
picture of the DOS would perhaps suffice to describe the 
evolution on going from Ti through Mn. The t2g mani- 
fold possesses a peaked region at the highest energy in all 
four compounds. In the Cr compound, this peaked region 
coincides with the Fermi energy. Interestingly, LiMn2 04 
seems to be a low-spin system, with the electronic config- 
uration i2g(T)j ^29^(1) J which is distinctly different from 
the manganese oxides perovskites that display colossal 
magnetoresistance jSSj, l39| . In the perovskites, there is 
a gap between filled t2g(T) states and the partially filled 
states which are mostly eg(l') with some small admixture 
of O p and ^2^(1). 

The compound ZnV204 which is isostructural (spinel) 
and nearly isoelectronic {d electron count of 2 on the B 
site) with LiCr204 is known to be a Mott insulator 



ordering antiferromagnetically at 45 K jij, associated 
with a tetragonal distortion of the spinel unit cell. Re- 
cently Reehuis et al. "4^ have performed a careful neu- 
tron diffraction study on this system. They find a struc- 
tural transition from cubic Fd3m to tetragonal /4i / amd 
at 51(1) K which breaks the magnetic frustration of the 
spinel structure and allows an antiferromagnetic ground 
state to be obtained (with Tn — 40(2) K). These authors 
have obtained precise nuclear and magnetic structures at 
60 K (cubic, paramagnetic) and 1.8 K (tetragonal, anti- 
ferromagnetic). We have performed LAPW calculations 
on the cubic phase, and LMTO calculations on using 
both the cubic as well as the low temperature tetragonal 
structure The LAPW structure optimization gave 

a cell parameter of 8.348(3) A to be compared with the 
60 K experimental value of 8.4028(4) A. The calculated 
internal x parameter of O was 0.259(4) A (experiment = 
0.2604(2)). The very close correspondence in this system 
provides confidence for the optimized LiCr204 structure. 
A magnetic moment of 1.84 fiB per V was obtained from 
LMTO calculations on the cubic phase (1.85/irmB from 



7 



CD 
O 



Q. 
W 

> 
CD 

tn 

CD 
■I—' 
CO 

-4—' 

-55- 
c/3 
O 
Q 



40 

30 

20 

10 


20 

10 



10 

20 



— metal d states 

— Op states 




1 

Cr 


'1 






1 


,1 







_8 -6 -4-2 2 
energy (eV) 



20 
10 


10 
20 





1 1 1 










J 




1 1 1 1 1 


l| ' 


! - 


8 -6 


-4 -2 


2 


4 



1 1 1 1 1 
Mn 1 


1 ' 








-8 -6 -4 -2 2 4 
energy (eV) 



FIG. 6: LMTO (partial) density of states for LiM2 04 (M= Ti, V, Cr, Mn). Black traces indicate M d states and grey traces 
indicate O p states. Both LAPW and LMTO calculations yield spin-polarized ground states for the compounds with M = V, 
Cr and Mn. 



LAPW), indicating some reduction of the moment from 
the spin-only value of 2 ^-q- The DOS for V d states and 
O p states for spin-polarized cubic ZnV2 04 is displayed 
in Fig. 7(a). Key differences in the DOS of ZnV204 and 
LiCr204 include reduced metal-oxygen covalency corre- 
sponding to fewer V d states in the region of O p states. 
In addition, the Fermi energy no longer lies on a sharp 
DOS peak, despite similarities in the shapes of the i2g(T) 
manifolds. Cubic ZnV204 is a poorer half-metal when 
compared with LiCr204 as already manifest in the re- 
duced magnetic moment. The d density of states of one 
of the V atoms in antiferromagnetic ZnV204 is displayed 
in Fig.[7|[b). The LMTO calculations were performed on 
the experimental structure explicitly considering four in- 
dependent V atoms forming a distorted tetrahedron with 
a short and a long V-V distance. The antiferromagnetic 
structure corresponds to spins on the tetrahedron being 
antiparallel when proximal (2.962 A) and parallel when 
far (2.963 A) in correspondence with experiment [43l |. 
We observe that the DOS remains metallic, indicating 
that correlation must play an important role in this com- 



pound, and that the one-electron description does not 
suffice. The changes in the DOS on going from the ferro- 
magnetic cubic structure to the antiferromagnetic tetrag- 
onal structure include a slight narrowing of t2g states and 
a reduction in the contribution from the DOS peak just 
below the E-p- 



E. Magnetism of LiCr204 

Spin polarized LMTO calculations suggested that the 
magnetic ground state for LiCr204 is more stable than 
the non-magnetic ground state obtained from a non-spin 
polarized calculation, by 0.4 eV per Cr atom. The con- 
verged magnetic moment on Cr was 2.49 /ib in calcula- 
tions that made use of the non-local Perdew-Wang ex- 
change correlation prescription. Since there are 2.5 d 
electrons per Cr"^-^+ in LiCr204, the value of the mag- 
netic moment corresponds to nearly 100% spin polariza- 
tion (an extensive discussion was given in section II B. 

Withing the usual von-Barth-Hedin LSDA the 
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FIG. 7: (a) Partial DOS for cubic ZnV204. Black traces 

indicate V d states and grey traces indicate O p states, (b) V 
d states for a single V atom in antiferromagnetic, tetragonal 
ZnV204. 



cell parameter a and the structural parameter 5 which is 
a measure of the trigonal distortion of CrOe octahedra as 
described in section II. Prom Table 1, we see that larger 
unit cells favor a larger magnetic moment and larger P. 
This is due to the narrowing of bands as the separation 
between atoms becomes larger. Large values of S seem 
to be contraindicated for magnetism and half-metallicity; 
the best half-metals correspond to nearly perfect CrOe 
octahedra. These rcsiilts were confirmed by the structure 
optimization procedure of the LAPW calculations. 
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TABLE 1: Magnetic properties of LiCr2 04 (LMTO/GGA) as 
a function of the crystal structure, a is the cell parameter in 
A, 6 is the deviation of the oxygen position from 0.25, M is 
the magnetization per Or in /ub , and P is the percentage spin 
polarization of conduction electrons. 



value was reduced to 2.42 /ib per Cr. We were interested 
in how changes in the crystal structure might affect the 
magnetism and half-metallicity in LiCr204. To this end, 
we have mapped in Table 1, the magnetization M per 
Cr, as well as the extent of half-metallicity at E-p defined 

N^{Ep)+N^{Ep)''^^^^'' 
as a function of the two structural parameters, the cubic 
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